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We derive the second-order dissipative relativistic hydrodynamic equations in a generic frame 
with a continuous parameter from the relativistic Boltzmann equation. We present explicitly the 
relaxation terms in the energy and particle frames. Our results show that the viscosities are frame- 
independent but the relaxation times are generically frame-dependent. We confirm that the dissi- 
pative part of the energy-momentum tensor in the particle frame satisfies JT*^ = obtained for the 
first-order equation before, in contrast to the Eckart choice ST'^" u„ = adopted as a matching 
condition in the literature. We emphasize that the new constraint 5T'^^ — can be compatible with 
, the phenomenological derivation of hydrodynamics based on the second law of thermodynamics. 

■ I. INTRODUCTION 

After the discovery that perfect hydrodynamics can be valid for describing the phenomenology of Relativistic Heavy 
. Ion Collider (RHIC) at Brookhaven National Laboratory [ll-Q , people are now interested in relativistic hydrodynamics 
Qh' for dissipative systems; see the recent excellent review articles [J, |5|. 

Q^, Recently, Tsumura, Kunihiro (the present authors) and Ohnishi (abbreviated as TKO) [6| derived generic covariant 
' hydrodynamic equations for a viscous fluid from the relativistic Boltzmann equation in a systematic manner with 
f~i , no heuristic arguments on the basis of the so-called renormalization group (RG) method [7Hll|. Although the hy- 
drodynamic equations they derived are the so-called first-order ones, the equations have remarkable aspects: The 
' generic equation derived by TKO can produce a relativistic dissipative hydrodynamic equation in any frame with an 
^ appropriate choice of a macroscopic flow vector (/z = 0, 1, 2, 3), which defines the coarse-grained space and time; 
^\ the resulting equation in the energy frame coincides with that of Landau and Lifshitz [l^] , while that in the particle 
frame is similar to, but slightly different from, the Eckart equation [13]. 

Let ST^'^ and SN'^ be the dissipative term of the symmetric energy-momentum tensor and the particle-number 
, vector, respectively. Owing to the ambiguity in the separation of the energy and the mass inherent in relativistic 
\^ • theories, one must choose the local rest frame (LRF) where the flow velocity u'^ with = 1 is defined: One 

[ of the typical frame is the energy (Landau) frame in which ST'^'^ A^p — with A''^ = g'^'^ — and g'^'^ = 

diag(-|-l, —1, —1, —1), i.e. there is no dissipative energy flow. On the other hand, another typical frame is the particle 
(Eckart) frame in which SN^^ A^^ — 0, i.e. there is no dissipative particle flow. Both in the energy and particle frames, 
the dissipative terms of the energy-momentum tensor and the particle-number vector are usually assumed to satisfy 
the constraints. 



X 

u^STf^^u^^O, (LI) 

and SN^^ = 0. These phenomenological ansatz have been eniployed as the matching conditions even in the 
subsequent "derivations" of the so-called second-order equations [14l4l7| : note that even in the Grad's moment method 
[isj , some ansatz are needed to 6T'^'^ and SN^^ as the matching conditions, for which different proposals exist [15, l^j. 

Here we emphasize that the matching conditions touch on the fundamental but not yet fully understood problem 
how to define the LRF in the relativistic fluid dynamics for a viscous system. The way how to deflne the LRF or 
equivalently to flx the matching condition is unsolved yet, and remains a nontrivial and fundamental problem in 
the field of nonequilibrium relativistic dynamics, although there have been no serious consideration on this difficult 
problem in the literature. Actually, we shall argue that these phenomenological ansatz, especially Ea. (ILl[) . can be 
false and actually is not compatible with the underlying kinetic equation. 

In fact, it is found that the TKO equation in the particle frame derived from the relativistic Boltzmann equation 
satisfies 

-^n 0' (1-2) 

but does not satisfy Eq. (II.l[) . One should here note that the derived condition (11.21) is identical to a matching condition 
postulated by Marie [1^ and advocated by Stewart [20[ in the derivation of the relativistic hydrodynamics from the 
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relativistic Boltzmann equation with use of the Grad's moment method. In their paper TKO proved that the 
Eckart constraint (jl.ip in the particle frame cannot be compatible with the underlying relativistic Boltzmann equation 
for the first-order hydrodynamic equation. In spite of the first-order one, the TKO equation in the particle frame is free 
from the pathological properties [2ll| in contrast to the original Eckart equation with which the thermal equilibrium 
becomes unstable for a small perturbation [22,]. 

One may naturally ask if the Eckart constraint p.ip should be replaced with (jL2]) even for the so-called second- 
order equation like Israel- Stewart one. And are any modifications needed to the constraints in the Landau frame? 
A purpose of this Letter (23j is to answer these questions both by phenomenological and microscopic analyses. We 
shall see that the Eckart constraint should be replaced with the new one even in the second-order equation, while 
no modification is necessary for the constraints in the energy frame. We shall derive the second-order dissipative 
relativistic hydrodynamic equations in a generic frame with a continuous parameter 9 from the relativistic Boltzmann 
equation. We shall derive the relaxation terms for a generic frame with the new constraint, and present explicitly 
those in the energy and particle frames. We shall show that the viscosities are frame-independent but the relaxation 
times are generically frame-dependent in accordance with the observation by Betz et al. |17| , although the constraint 
to ST^"' is quite different. 



II. A GENERAL PHENOMENOLOGICAL DERIVATION OF RELATIVISTIC DISSIPATIVE 
HYDRODYNAMIC EQUATIONS; EXISTENCE OF POSSIBLE EXTRA TERMS IN THE DISSIPATIVE 

TERMS 

Let T'^'^ and N'^ be the symmetric energy-momentum tensor and the particle-number vector of the system we 
consider, respectively. The total number of independent variables is fourteen, and the dynamical evolution of these 
variables are governed by the respective balance equations; 

a^T^" - 0, (ii.i) 

^^N^' = 0. (II.2) 

With use of an arbitrary four vector with = 1^ T^" and can be cast into the tensor-decomposed forms, 

T^"" = {e + Se)u^'u'' -{p + Sp)A^"' + q^'u'' + q''u^' + TT^''', (II.3) 
= + (II.4) 

respectively. Here, e -I- Se, p -\- dp, and n + Sn are the internal energy, pressure, and particle-number density in the 
dissipative system; e + Se = Tab u°' u^, p + Sp = —1/3 Tab A"'', and n + Sn = Na u°', with e = e{T, n), p = p{T, fj,), 
and n = n{T, fi) being the corresponding quantities in the local equilibrium state characterized by the temperature 
T and the chemical potential /i. Note that we have made it explicit by Se, Sp, and Sn that the dissipations may cause 
corrections to all these quantities, although only the correction to the pressure has been considered in the literature; 
Sp is identified with the bulk pressure 11. We emphasize that there is no persuading reasoning that only the pressure 
acquires corrections due to the dissipative process. The dissipative parts of the energy- momentum tensor and particle- 
number vector are identified as ST^'^ — Se u'^ + Sp A'^" + u'^ + (f -I- -n^'^ and SN^ = Snu^ + , respectively. 
The energy flow relative to is denoted by q^ , is the flow of particle number relative to , and flnally tt^" is the 
shear stress tensor; q^ = Tab u° A''^, i/^^ = Na A"^, and tt'^'^ = Tab A"''^''. Here the space-like, symmetric and traceless 
tensor At"'P'^ = l/2{Af'PA'"' + Af"" A^p - 2/3 A^" AP") is introduced. One can easily confirm that qP Uf, = 0, 
fP Ufj^ = 0, ttP"^ = tt'^p, and u^-kP'^ = ttP^ = 0. This implies that the total number of independent components of 
qP, i/P, and ttP'^ is eleven. Since TP'^ and A^'^ have the fourteen components in total, Se, Sp, and Sn have only one 
independent component other than T and ^. We take Sp = H as the independent component as a natural choice, then 
Se and Sn can be expressed as Se = feU and Sn — /„ H, where /e and /„ are functions of T and fj,; /e = fe{T, /i) 
and /„ = fn{T, jj). Here we have assumed that the dissipative order of Se and Sn are the same as that of Sp at most. 
We remark that although /e and /„ may take finite values generically, the functional forms of /e and /„ cannot be 
determined by the phenomenological theory, as those of e, p, and n can not, either. All the previous analyses assumed 
fe = fn = 0, which has not been recognized so far. 

Now we shall show that the just usual phenomenological derivation of the hydrodynamic equations in which the 
second law of thermodynamics is utilized allows the existence of Se and Sn, i.e., finite values of /e and /„, in the 
relativistic dissipative hydrodynamic equations. It is found that the essential point of the proof is the same for the 
first- and second-order equations where the entropy current Sp is at most linear and bilinear with respect to H, qP, 
fP, and TT^j^, respectively, although the resulting mathematical expressions are much more complicated in the second- 
order one [2J| . Thus we here take the first-order equation, for the sake of simplicity. The second-order equations with 
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finite /e and /„ will be derived microscopically later in this article. So the entropy current is given by 

T S^' =pu'' + u^T'"' - ^lN^'. 

The second law of thermodynamics reads d^S^ > 0. 
The divergence of S"^ is found to take the form 



(II.5) 



n 



feD 



1 



1 



T T 



7^ Du,, 



(II. 



where D = u"' da and = A'^''9a. Here, we have used the conservation laws, Eq.'s (III.ll) and pL2p . and the first 
law of thermodynamics, D{p/T) + eD{l/T) - nD{fi/T) = 0. 

The frames define the flow velocity u'^ of the fluid: The flow velocity in the particle frame and the energy frame 
are defined by setting u'' = Nf^/^N" and m'' = T^" Ua/^/T"'' UbT^^c respectively [l^. By these s ettings, 
a closed system of the relativistic dissipative hydrodynamic equations is obtained. Note that = N'-'- / \/N^ N^, 
(mM = T^" Ua/\/T''^UbT^c u") is equivalent to v^' = Q (g^ = 0). 

In the particle frame where = 0, Eq. pi.6p is reduced to 



a^^" = n 



— DUn + Vu — 



(11.7) 



It is found that the following constitutive equations, 



n 



-XT'" 



T 



T 



(II.8) 

(11.9) 
(11.10) 



guarantees the second law of thermodynamics, d^S^^ > 0, with C, A, and t] being the bulk viscosity, heat conductivity, 
and shear viscosity, respectively. This is because the divergence 9^5^ now becomes positive semi-definite; 



(T 



2r]T 



> 0. 



(11.11) 



Thus we realize that there is nothing wrong with the resultant relativistic dissipative hydrodynamic equations with 
finite /e and /„, or equivalently finite 6e and 6n, which is compatible with the second law of thermodynamics. Eq.'s 
(III.8l) - (jII.10|) with a restricted condition fe = fn — are identical to the constitutive equations proposed by Eckart 
that are commonly used. 

In the energy frame where = 0, we can obtain the constitutive equations in the same way as the particle-frame 
case with /e and /„ being kept finite. The resultant equations are given by Eq.'s (jll.8[) . (jII.10[) . and 



T' 



(11.12) 



with h = (e + p)/nT being the enthalpy. It is noted that these equations are reduced to the constitutive equations 
by Landau if we can set fe — fn ~ 0- 

By applying the above argument to the entropy current at most bilinear with respect to 11, q^, v^^, and tt^'', we 
can obtain the relaxation equations with /e and /„ being finite, which make up the so-called second-order relativistic 
dissipative hydrodynamic equations together with the conservation laws in Eq.'s pi.l[) and (III. 21) [24]. 

Now the dissipative part of the energy-momentum tensor satisfies ST'^'^ = Se = feU and 5T'^^ = Se — S 6p = 
( fe — 3) n. As emphasized before, the values of fe and /„ can be determined only from a microscopic theory. The 
phenomenological theory cannot proceed further because no such logic to determine them is implemented in the 
theory. In the following section, we shall show that the microscopic theory gives /e = 3 together with /„ = in the 
particle frame while fe — fn — in the energy frame, and hence ST^^ = but ST'^'^ u^, = 3 11 ^ in the particle 
frame. This fact tells us that the usual constraint employed for the particle frame must be abandoned, and all the 
analyses based on this constraint should be redone. 
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III. DERIVATION OF SECOND-ORDER EQUATIONS AS LONG WAVELENGTH AND LOW 
FREQUENCY LIMIT OF RELATIVISTIC BOLTZMANN EQUATION 

The argument so far is in the stage of thermodynamics where the argument is robust but the parameters such as 
fe and /„ as well as the equations of state e, p and n appearing in the theory remain undetermined. The problem 
which we encounter is how to reduce a dynamical equation to a slower one described with fewer dynamical variables. 
For this purpose, we will investigate the infrared limit of the relativistic Boltzmann equation with use of a powerful 
reduction method, the "RG method" 0-[ni- 

The RG method is a systematic reduction theory of the dynamics leading to the coarse-graining of temporal and 
spatial scales. The full presentation of the reduction of the relativistic Boltzmann equation to the second-order 
hydrodynamic equation is technical and involved. So we here only present main results with key several equations, 
leaving the detailed account to another publication [24], although the derivation of a wide class of the first-order 
equations is presented in Ref.f3|. 

We start with the simple relativistic Boltamann equation, 

p^^d,fp{x)^c[fUx), (iii.i) 

where fp{x) denotes the one-particle distribution function defined in the phase space (a;'^ , p'^) with p^^ being the 
four momentum of the on-shell particle. The right-hand side of Eq. pil.ip is the collision integral, C[f]p{x) = 
m ^ T,p2 ^ Pi b2 , Ps) (/p2(a^)/p3 (2;)- /p (2^) /pi (a;)), where a;(p, Pi |p2 , Pa) denotes the transition 

probability owing to the microscopic two-particle interaction. 

We are interested in the hydrodynamical regime where the time- and space-dependence of the physical quantities 
are small. In another word, the time and space entering the hydrodynamic equation are the ones coarse-grained from 
those in the kinetic equation. Thus we are lead to introduce a macroscopic Lorentz vector, a^{x) which specifies the 
covariant coordinate system and we call the macroscopic-frame vector. With use of (x) , we define the macroscopic 
covariant coordinate system (r , cr'') as dr = a^{x) dx^ and dc'^ = {g'^'^ — al^{x)ap{x)/ap{x)) dx^, = A^^lx) dx^. 
We note that the small quantity e has been introduced to tag that the space derivatives are small for the system we 
are interested in. e may be identified with the ratio of the average particle distance over the mean free path, i.e., the 
Knudsen number. 

In this coordinate system, Ea. (|III.ip can be cast into 

|-/p(t , a) = ] C[/]p(T ,a)-s 1 p ■ V/p(r , a), (III.2) 

where a^{T , a) = al^{x), A^''(r,cr) = A^^^x), and fp{T,a) = fp{x). Since e appears in front of = 
A^"^(r , a) Ea. pil.2l) has a form to which the perturbative expansion with respect to e can be applied. In the per- 
turbative expansion, we shall take the coordinate system where a^(T , a) has no r dependence, i.e., a^(T , a) — a^{a). 

The zeroth-order approximate solution we construct is a stationary solution, which is identical to a local equilibrium 
distribution function given by the Juetner function /^^ = {2tt)~^ exp[(/i — p • u)/T]. Note that this solution contains 
five would-be integration constants, T, /x, and with u^^ — 1, which can be identified with the temperature, the 
chemical potential, and the fluid velocity, respectively. 

The collision integral is expanded around the zeroth-order solution and is reduced to the linear operator Apq = (p • 
"p)~^ WT^if^^p- Furthermore, it is found to be convenient to convert Apg to Lpq = fp^~^ Apq f^^ = [/°^~^ Af'^'^]pq, 
with the diagonal matrix = fp'^Spq. We also define the inner product between arbitrary vectors ipp and ipp by 

(^, ^(?'-°p)/p>pV'p- (ni.3) 

p ^ 

With this inner product, we can define a normed linear space. 

Now the first-order solution is given in terms of the five zero modes of L, (fQp = (1, p^). The corresponding variables 
are just T, fi, and with u^m^ = 1. The zero modes span a linear space Pq, which is an invariant manifold for the 
asymptotic dynamics of the relativistic Boltamann equation in the terminology in the dynamical systems 0, [25| . 

Then the second-order solution is given by incorporating the next slow modes, which span a linear space Pi. We 
naturally require Pi is orthogonal to Po, that is, Po -L Pi. We find that Pi is expanded by the bilinear forms 
of momenta; (p^p = [Qo'f^'^]p, where (p!^'^ = p^ p'^ , and Qq is the projection to complement to Po- By definition, 
( ip^'^ I ) = is satisfied. Note that the dimension of (p^p is nine, which correspond to the number of the new 
would-be integration constants, 11, J'^ with J'' — 0, and tt'^'^ with tt^" = tt"^^ and tt^'' — tt^^ — 0. 
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A generic choice of the macroscopic frame vector is = {{p ■ u) cos 6 + m sin0)/(p ■ u)u^, where is a parameter 
defining the frame. For example, 6 — Q (9 = tt /2) gives the energy (particle) frame. 

The resultant generic relaxation equations of the second-order hydrodynamic equation with 9 being kept are 

n = Xn-mDIi-tnjV\]a+Xnii'R + X^jJa+X^,TTab, (III.4) 
= - Tj A^"^ DJa - l.m V^n - Ij^ A^''"- Va^bc + ^^n n + J a + ^ab, (ULS) 
= - A^"^"" D^ab - 4 J A^'"^^ VaJb + ^^n n + X^J J„ + X^^' ^ab- (111.6) 

Here, Xn, Xj, and X!^^ are the thermodynamic forces; their simple forms retaining only Xu, Xj, and X^'^ are the 
usual constitutive equations. The relaxation equations of 11, J^, and tt'"^ are characterized by the relaxation times 
Tn, Tj, and Ttt, while €njj ^jiIj ^j-n-i and f^rj mean the relaxation lengths. The correction to the thermodynamic forces 
Xu, X^j, and X^- are given by Xun, X^j, X^l X^}^, Xfj, X^f , X^?^, and X^.^^ ^ 

The continuity equations of the second-order equation in the energy frame is found to be given by setting = as 
in the first-order case @ and read d^T'^'^ = and df^N^^ = 0, where 

Tf"" = ew^'u"^ - (p + n) A'"' +71^"^, (III.7) 
N^' = nu'' + J>'. (III.8) 

The thermodynamic forces are Xu = -C Vua, ^ Xh-'^ W^'{n/T), and X^" = 2 ry A'"'"'' VaUb, which clearly show 
that fn=0 and /e = as was anticipated. 

The energy- momentum tensor and particle-number vector in the particle frame with 6 = tt/2 read 

T^'^ = (e + 3 n) u'' - (p -t- n) A^"^ -I- J"^ -t- w'' -I- tt^'', (III.9) 
iV^ = nu'^, (III.IO) 

and Xn -C(3 7- 4)-^ (v°m„ -3TDT-^), X^ = -AT^ (V^T^i -I- T^^ Du^), and = 2 r/ A^'^'^'' VaUb, where 
7 = 1 + (z^ — /i^ -I- 5 /i— 1)^^ is the ratio of the specific heats. Thus we find that /e = 3 with /„ = 0, as we announced. 

Although we have obtained the relaxation equations for the dissipative forces 11, J'^, and tt'"^ for arbitrary 9 p3 |. 
we shall only write down them for two typical frames, i.e., the energy {9 = 0) and the particle {9 = tt/2) frames for 
the sake of the space. 
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(A) In the energy frame = 0): 
n = -CV"u„-Tn£'n-^njW„ 



^^n. 4°)v«^-.«i?.« + g4(|HA- 



(III.ll) 



n 



^r. A- 



- i I A-"" Ve^ - «W ^-c) + ^ A'^, A^<^«/) A,/"! 
tt'*'' = 2ryA''^«''VaUb-4jA'*''«^VaJ6-r,A''^»''£)7r„6 

-^4n|-«.nA'^'''^''VaW6|n 

- \ i.j I A-- - I^.,) + g A-,, A--) A/ 1 J„ 



(III. 12) 



(III. 13) 



where w'"' = {V^'u" - \7''u'')/2 is the vorticity. 
(B) In the particle frame {6 = 7r/2): 



n = -C(37-4)-2(v'^u„-3TI?i) -7ii£)n-fnjWa 



1 I ™ C(37-4)-2T^ 

7^ rn <^ Kn + -^^^ ^ da 

i Tn 



.C(37-4)-2T 



1 



(")v7a^ .,(1) ^..g , C(37-4)-^r 



nj - 



A^^ A/ >> J„ 



C(37-4)-2r 



1 1 



(III. 14) 



- \ Am - Du^ + ^ A^„ A 



/n A"M J,n 



- \ ij^ (a-- (4°) - 4^i ^-c) + ^ A-/) A,/4 



(III. 15) 
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FIG. 1: The 9 dependence of Tn and tj at m/T — 0.5 and fi/T — 0.0. We normalized the relaxation times by the corresponding 
transport coefficients. The energy and particle frames correspond to 6 = and 7r/2. 



and Eq. (IIII.13j) . Note that the effective bulk viscosity Coff = C (3 7 - 4)"^ [lH appears in Eq. (|III.14p . 

Here we have introduced the new coefficients, Iutt, i-n-n, ^n, '^nj; '^njj i^Utt, 'Sjni '*jn: \ '^/'ii '^/^i Kirn, 



At^^j, K^j, ti^'', and K^-^\ which are complicated functions of T and ^ [24 1 



We have found that the relaxation times are frame dependent except for r^r while the transport coefficients such 
as the viscosities and the thermal conductivity are frame independent. For a demonstration of the frame-dependence 
of the relaxation times, we show in FfGll]the frame (9) dependence of rn and tj, which tends to increase when the 
frame changes from the energy to particle frame. 



IV. BRIEF SUMMARY 



In summary, we have derived the second-order dissipative relativistic hydrodynamic equations in a generic frame 
with a continuous parameter 9; the generic frame is reduced to the energy and particle frame with the parameter 
choice 9 — and 7r/2, respectively. A notable point of our result is that the dissipative part of the symmetric 
energy- momentum tensor ST^''' in the particle frame satisfies the equality (5r^ = 0, in contrast to the usual choice 
Uf^ST^^Uiy — 0, while STf^" of our derived equation in the energy frame satisfies the usual constraint u^j^dT^" Ui, = 
0. We emphasize that this novel equality in the particle frame is a consequence of the derivation based on the 
renormalization-group method, a powerful method for the reduction of dynamical systems. We note that the same 
constraints were also derived for the first-order dissipative relativistic hydrodynamic equation [2l| . We have also 
shown that the phenomenological derivation based on the second law of thermodynamics allows that ST^^'^ can 
be proportional to the bulk pressure 11 and non-vanishing in the particle frame. Indeed, our microscopic derivation 
shows that ST'^" = 3 11. We have presented the relaxation equations in the energy and particle frames, expli citly 
as typical examples, although we have obtained the microscopic expressions for them in a more generic frame [24| . 
We have shown that the viscosities are frame-independent but the relaxation times are generically frame-dependent, 
as depicted in FIG[T] The detailed derivation of the equations and discussions on the phenomenological consequences 
of the hydrodynamic equations thus obtained will be discussed in forthcoming papers [2^ . 



acknowledgment 

We thank Tetsu Hirano for his interest in this work and discussions. T.K. thanks Dirk Rischke for his interest in 
our work. This work was partially supported by a Grant-in- Aid for Scientific Research by the Ministry of Education, 
Culture, Sports, Science and Technology (MEXT) of Japan (No. 20540265), by Yukawa International Program for 
Quark-Hadron Sciences, and by the Grant-in-Aid for the global COE program "The Next Generation of Physics, 



8 



Spun from Universality and Emergence" from MEXT. 



[1] See review articles, P.Huovinen, in "Quark Gluon Plasma 3", ed. R.C.Hwa and X.N.Wang, (World Scientific, Singapore), 

p.600; P.F.Kolb and U.W.Heinz, ibid, p. 634. 
[2] T. Hirano and K. Tsuda, Phys. Rev. C 66, 054905 (2002); D. Teaney, Phys. Rev. C 68, 034913 (2003). 
[3] M.Gyulassy and L.McLerran, Nucl.Phys.A 750, 30 (2005). 
[4] T. Hirano, N. van der Kolk and A. Bilandzic, arXiv:0808.2684 [nucl-th]. 
[5] P. Romatschke, arXiv:0902.3663 [hep-ph]. 

[6] K. Tsumura, T. Kunihiro and K. Ohnishi, Phys. Lett. B 646, 134 (2007). 

[7] L. Y. Chen, N. Goldenfeld and Y. Oono, Phys. Rev. Lett. 73, 1311 (1994); Phys. Rev. E 54, 376 (1996). 

[8] T. Kunihiro, Prog. Theor. Phys. 94, 503 (1995); 95, 835 (1996) (E); Jpn. J. Ind. Appl. Math. 14, 51 (1997); Prog. Theor. 
Phys. 97, 179 (1997). 

[9] S.-I. Ei, K. Fujii, and T. Kunihiro, Ann. Phys. 280, 236 (2000). 
[10] Y. Hatta and T. Kunihiro, Ann. Phys. 298, 24 (2002). 
[11] T. Kunihiro and K. Tsumura, J. Phys. A 39, 8089 (2006). 

[12] L. D. Landau and E. M. Lifshitz, Fluid Mechanics, (Pergamon Press, London, 1959). 

[13] C. Eckart, Phys. Rev. 58, 919 (1940). 

[14] W. Israel, Ann. Phys. 100, 310 (1976). 

[15] W. Israel and J. M. Stewart, Ann. Phys. 118, 341 (1979). 

[16] S. R. de Groot, W. A. van Leeuwen and Ch. G. van Weert, Relativistic Kinetic Theory, (Elsevier North-Holland, 1980). 
[17] B. Betz, D. Henkel and D. H. Rischke, arXiv:0812.1440 [nucl-th]. 
[18] H. Grad, Comm. Pure Appl. Math. 2, 331 (1949). 

[19] C. Marie, Annales de I'institut Henri Poincare (A) Physique theorique, 10, 67 (1969); 127 (1969). Marie's theory is reviewed 
in [13, [11. 

[20] J. M. Stewart, N on- Equilibrium Relativistic Kinetic Theory (Lecture Notes in Physics No. 10; Springer, Berlin, 1971). 
[21] K. Tsumura and T. Kunihiro, Phys. Lett. B 668, 425 (2008). 
[22] W. A. Hiscock and L. Lindblom, Phys. Rev. D 31, 725 (1985). 

[23] K. Tsumura and T. Kunihiro; talk presented at JPS meeting at Yamagata University, September 20, 2008; T. Kunihiro, 

Y. Minami and K. Tsumura; talk presented at Quark Matter 2009, Knoxville, March, 2009. 
[24] K. Tsumura and T. Kunihiro, in preparation. 

[25] See for example, J. Guckenheimer and P. Holmes, "Nonlinear Oscillators, Dynamical Systems, and Bifurcations of Vector 
Fields" Springer- Verlag, 1983. 

[26] C. Cercignani and G. M. Kremer, The relativistic Boltzmann equation: theory and applications, (Birkhauser, 2002); Progress 
in mathematical physics v. 22. 



